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I. INTRODUCTION 



Quantum walks are the simplest formal analogues of classical random walks. They have 
been introduced in their discrete-time version by^ and^^ and the continuous-time version first 
appeared in^. Quantum walks are important in quantum information and quantum com- 
puting, where several algorithms are based on them; examples include Ambainis's algorithm 
for element distinctness and an algorithm for the triangle finding problem^-. Quantum 
walks are also useful in several physical or biological contexts. They are naturally systems 
of choice to explore fundamental issues in quantum physics, including decoherencei^'2^, but 
they also prove useful in modeling transport in solids^ii"^^ disordered media^i^i^ and even 
complexes of algae^^'i^. Quantum walks have been realized experimentally, for example as 
transport of trapped ions^SiM^ of photons in wave guide lattices^^ or optical networks^S and 
of atoms in optical lattices^. Quantum walk experiments of two photons^^ have recently 
been performed, with the possibility of simulating Bose or Ferni statistical^. Cavity QED 
quantum walks have also been proposed^I. 

The continuous limit of discrete time quantum walks with constant coefficients has al- 
ready been addressed by several author3^>^'^>22i>^, with very different and apparently not 
reconcilable approaches. This article revisits this problem and obtains new, definitive con- 
clusions for one-dimensional discrete-time quantum walks, which are entirely determined by 
three possibly time- and space-dependent Euler angles^. Our main results are the following: 
(i) generically, the object which can admit a continuous limit is not a given discrete-time 
walk, but rather a 1-jet of discrete-time walks, which regroups all walks which share a com- 
mon expansion in the time-and length steps as these tend to zero (ii) only a few 1-jets admit 
a continuous limit (iii) when it exists, the continuous limit can be described by a 2D Dirac 
-like equation or, alternately, two independent Klein-Gordon equations (iv) the Dirac-like 
equation derives from a variational principle which exhibits an interesting symmetry. 

II. FUNDAMENTALS 

A. Discrete-time quantum walks 

We consider quantum walks defined over discrete time and discrete one dimensional space, 
driven by time- and space-dependent quantum coins acting on a two-dimensional Hilbert 
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space "H. The walks are defined by the following finite difference equations, valid for all 

(j, m) G N X Z: 



where 



B 



1 ^j,mi C 



.j,m.J 



m—l 



e*^ COS 9 e*'' sin 9 
-e~*^ sin 9 e~*^ cos 9 



(1) 



(2) 



is a SU(2) operator defined in terms of its three Euler angles 9,C,,C- The index j labels 
instants and the index m labels spatial points. The wave function has two components 
ip~ and which code for the probability amplitudes of the particle jumping towards the left 
or towards the right. The total probability tt^ = J2m (I ^j~m P + I ^/m P) independent of 
j i.e. conserved by the walk. The set of angles {9j^m,^j,m,Cj,m, (j,^) G N x Z} defines the 
walk and is, at this stage, arbitrary. 

To investigate the continuous limit, we first introduce a time step At and a space step Ax. 
We then introduce, for any quantity a appearing in ([T]), a function a defined on : 



that the number aj^m is the value taken by a at the space-time point [tj = jAt, Xr, 
Equation ([1]) then reads: 



M such 
niAx). 



^ {tj + At, Xm) 
ij+{tj+At,Xm) 



B (9(tj^ ^^m)) ^(j^j: •^m): C(^i! -^m)) 



1p {tj,Xm + Ax) 
ij^(tj,Xm - Ax) 



(3) 



where the tilde's have been dropped on all functions to simplify the notation. We now 
suppose, that all functions can be chosen at least in both space and time variables for all 
sufficiently small values of At and Ax. The formal continuous limit is defined as the couple 
of differential equations obtained from (eq:defwalk) by letting both At and Ax tend to zero. 

B. Scaling for the continuous limit 



Let us now introduce a time-scale r, a length-scale A, an infinitesimal e and write 

At = re 
Ax = Xe\ 



(4) 



where 6 > traces the fact that At and Ax may tend to zero differently. 
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For the continuous limit to exist, at least formally, the operator B{6, ^, () defined by the 
functions 6, C, and ( must also tend to unity as e tends to zero. This is so because the 
two column vectors of the left-hand side and on the right-hand side of (E]) both tend to 
\E'(tj, Xm) when At and Ax tend to zero. This trivial remark implies that the two functions 
6 and ^ must actually depend on e and tend to or tt as e tends to zero. In other words, it 
generally does not make sense to consider the continuous limit of a given walk, defined by 
e-independent angles, but rather the limit of a family of walks indexed by e, whose defining 
angles depend on e and tend to zero with this infinitesimal. Since we are only interested in 
the behavior of the family as e tends to zero, we only need to consider families which are 
jets and we thus write: 

9{t,x) =pTr + e{t,x)e" 

^{t,x)=p7i + ^{t,x)e^, (5) 

where p = or 1 and a, P > 0. Note that the third Euler angle ( does not have to tend to 
zero with e because of the sin^ factor in front of e*'' in ([2]). Note also that the Hadamard 
walk cannot have a continuous limit. This point is further discussed in Section 5 below. 

The continuous limit can then be investigated by Taylor expanding ilj^{t,x^Ax), ^/'^(t ± 
At, x), cos^, sin^ and e*^ in powers of e. The possible scalings obeyed by the continuous 
limit are found by examining the lowest order contributions. One has: 

ip^{t ± At, x) = ip^{t, x) + 0(e), 

V'^(t, X T Ax) = V^(t, x) + O(e^), 

e'^cose = l + 0(e^) + 0(e2"). 



(6) 

(7) 
(8) 



e'^ sine = 0(e' 



+ 0(n + Oie") + 0(e 



(9) 



(10) 



and 

Equation ([3]) then leads to: 

+ 0(e) = 

ij+{tj,Xm)j ^i)+{tj,Xm) 

Zeroth order contributions cancel out as expected and the remaining terms must balance 
each other. The richest and most interesting case corresponds toa = P = 6= l because 
all contributions to ( ITOj) are then of equal importance. This is the scaling that will be 
investigated in the remaining of this article. 



III. EQUATIONS OF MOTION 



The natural space-time coordinates to investigate this scahng are the so-called null coor- 
dinates u~ and u'^, defined in terms of x and t by 



1ft X 



, 1ft X 



(11) 



Partial derivatives with respect to these null coordinates read: 



d- = a,- = rdt - Xd^ 



d+ = du+ = rdt + Xd^, 



and 



= r'du - X^d,, = X^n, 
where □ is the usual d'Alembert operator, defined by 

n — -^dtt — dxx, 



(12) 



(13) 



(14) 



with c — X/t. 

The expansion of the discrete equations around e = then leads to the following equations 
of motion for ■0'" and ijj'^ : 



and 



d_ij;- = {i^^jj- - (-l)f+^^e+^4+) 



(15) 



(16) 



These coupled first-order equations are best transcribed in operator form i.e. as 7?^' = 0, 
with the operator V acting on the two-component wave function ^ given by: 



V = T-d- + r+a+ - iasi + ia2{-lY^^9 cos C + iai {-1)^+^9 sin C, 



(17) 



where the cr's are the three Pauli matrices: 



1 

1 



(18) 
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(19) 

(20) 

and 

= ^ (1 ± ^3) . (21) 

The operator V defines the formal continuous limit of the quantum walk. Equations (fT5|) 
and (fT6l) are sometimes called the Dirac form of the continuous quantum walk dynamics. 

These first-order equations imply that each component of \E' obeys, at points where 9 
does not vanish, an uncoupled Klein-Gordon equation. For example, the equation obeyed 
by ip' can be obtained by using (fT5|) to express ip^ in terms of ip~ and d-ip~ , and by then 
replacing by this expression in f|T6|) . This leads to 



0"2 



-i 

1 



(^3 



1 

-1 



X'^Dip- = [d+ {\ne + <) - d.ip- + i^d+ip- 

- [9^ + f + i^d+ (In 9 + iC)- id+^ ij-) (22) 

and 

X'^nilj+ = -ild.tl)^ + [(9_ (In^ - zC) + ii\ d+^p^ 

- + f + ild^ (- In ^ + iC) + id^i\ 7/-+.) (23) 

These equations ressemble Klein-Gordon equations, but contain terms which violate time- 
reversibility. The discussion of these equations is postponed till section 4.2. 

The natural initial conditions for the first order equations (|T5l) and ( JT6l) are simply the 
values of and ip^ at time t = i.e. on the submanifold u~ = —u'^. The natural initial 
conditions for (l22l) are the values of ip' and dtip' on the sub-manifold t = 0. By equation 
( !T5|) . giving oneself these initial conditions is equivalent to giving oneself both ^~ and 
at time t = 0. The same reasoning applies to the initial conditions of (123]) . 
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IV. VARIATIONAL PRINCIPLES AND SYMMETRIES 



A. First order form of the equations 



It is straightforward to check that the Lagrangian density: 



(24) 



leads to the first-order equations of motion (fT5 | [T6|) . This density and, hence the dynamics, 
is not invariant under the (global or local) action of SU{2). If it were, one would have 
VDU^^ = V i.e. [DjU] = for all elements U of SU(2). The generator o"3 does commute 
with P, but (Ti and a2 do not. This reflects the fact that cti and o"2 interchange ip' and ip^ 
while these two components of \E' play fundametally different roles, one propagating towards 
the left while the other propagates towards the right. This in turn is reflected by the fact 

that r+ ^ r-. 

The first-order dynamics can nevertheless be linked to interesting locally invariant oper- 
ators. Let us consider the following family of differential operators: 



where /x G {— ,+} is the space-time index, j G {1,2,3} and the S's are arbitrary real, 
possibly time- and position dependent quantities; the T matrices have been defined above 
in on]). 

It is obvious that the operator V characterizing the continuous dynamics of the quantum 
walks considered in this article belongs to this family. Let us now consider, for any j G 
{1, 2, 3} and arbitrary time- and space-angle a, the following transformation of T>{B): 




(25) 



V{B)j^o.{^) = exp {-iaaj) [I?(-B)(cosa^) + iajV{B) {sin a"^)] . 



(26) 



Replacing V{B) by its definition ( l25l) and using the fact that = 1, one finds: 



P(S),,,(^) = V{B){^)+za,{m,a)^. 



(27) 



In other words: 



V{B) 



(28) 



with 



B%a) = B'; + 5%a. 



(29) 
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The family of operators I){B) thus displays local invariance with respect to the transforma- 
tion fl26l) and the operator fixing the continuous quantum walk dynamics is one member of 
this locally invariant family. Let us stress that the operators V{B) do not commute with 
(Ti and (72"^ and that the invariance considered here is not an invariance of the continuous 
quantum walk under the action of SU{2). If it were, the aj in the bracket on the right-hand 
side of equation (!26l) would appear on the right of T>{B), not on the left. The invariance 
considered here is a formal invariance of an operator family, not an invariance of the equa- 
tions of motions. Note however that, for j = 3, the above transformation coincides with a 
standard gauge transformation along the ^3 generator of SU{2), because commutes with 
allV{B). 

Let us remark that a generic V{B) operator does not conserve the total probability dx[\ 
ilj~{t,x) p -|- I ilj~{t,x) p] for all ip' and If one imposes total probability conservation 
as an identity valid for all values of ip~ and ■?/'"'", one obtains that V{B) necessarily takes the 
form: 

V{B) =T^'df,-ia3{T^'Q+ia2{-iy^+^ecosC + i(ri{-iy+^esmC (30) 

where p is an arbitrary integer and 6 and ( are arbitrary possibly time-and position 

dependent real quantities. This coincides with the Dirac form of the continuous quantum 
walk, except for the fact that C,- is not necessarily identical to This can be remedied by 
exploiting the local U{1) invariance of the problem i.e. by changing the phase of \I'. The 
operator P given by equation (JTTll is thus the most general V{B) operator which preserves 
the total probability as an identity of the dynamics. By ( 129|) and (fT7|) . this property is 
preserved under a transformation f l26l) if and only if d^a = d^a i.e. d^a = 0. 

Let us end this section by writing out two useful forms of the Lagrangian density Ld- 
The first form serves as a reminder of the gauge invariance described above. It reads: 

Lcf^,^'^] = ^T^V;,^ (31) 

with 

= d^-ta3^+ta2{-iy^^ecosC + iai{-iy^^esmC (32) 

This can be derived directly from flMj) by noting that F" + r+ = L 

The Lagrangian density L^, can also be brought to a form which closely resembles the 
Dirac Lagrangian density. This is accomplished by singling out the C,ct3 term, which then 
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appears as a mass term (remember this term is related to the U{1) invariance of Let us 
define 

= -o-sr^ (33) 

and 

= d^ + iai{-iy+^esmC + ia2{-lY^^ecosC. (34) 

One can then write: 

Lz3[^,^^] = ^t(27^D^-C)^, (35) 

which does look like the standard Lagrangian density for a Dirac field. Of course, \1/ here 
belongs to a two -dimensional Hilbert space, and is not a Dirac spinor. 



B. Second order form of the equations 

Let us now investigate if each of the equations (1221) and (1231) can be obtained from a 
variational principle. As already noted, the second order equations fl^ and fl2^ are not 
time- reversible. This irreversibility is made manifest in (122]) by the {d^\n6)d-ip~ term and 
in (!23|) by the {d- In 6) d+il!~^ term. We will only present autonomous variational princi- 
ples for equations fl22l) and fl23l) in the special cases where these irreversible terms vanish. 
We nevertheless remark that a standard method for obtaining a variational principle for 
an irreversible dynamics is to enlarge the configuration space of the system and actually 
obtain a variational principle which delivers, not only the desired irreversible dynamics, but 
also another, coupled dynamics obeyed by some extra degrees of freedom. For the continu- 
ous quantum walks dynamics considered here, there are indeed no autonomous variational 
principles which deliver separately the Klein-Gordon equations obeyed by ip~ and when 
disspative terms do not vanish, but there is always a single variational principle which de- 
livers the coupled Dirac-like equations obeyd by the \E' = [tp' 

Let us start with equation fl22|) and suppose that disspative terms vanish i.e. that d^6 = 0. 
The only reasonable candidate to build a variational principle for ( l22l) is the standard Klein- 
Gordon Lagrangian density 

Lkg[^-\ = 9'''{y-r)^VZr - Ml I |2 (36) 
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where 

1 



(37) 



1 

are the metric components in the null coordinates {u^,u'^) and the covariant derivative 
= — iA^ couples tp^ to a real potential A~ . Note also that the mass M_ appearing 
in (1361) must be real for the Hamilton equations obeyed by ip~ and ip'"^ to be equivalent, as 
they should. The physical reason for this restriction on the mass M_ is that a mass with 
imaginary part would generate an imaginary part in the frequency (energy) or wave- vector 
(impulse) of the particle, and this would signal irreversibility. 

Identifying the equation of motion derived from (!36|) with equation (122|) leads to: 

M_ = \e\ (38) 

provided 

A2nC = (9_ + 9+)e = r9,e. (39) 

The only non possibly vanishing components of the electromagnetic tensor F (i.e. the 
curvature of the connection A) are F_+ = A straightforward calculation shows that 

these components actually vanish if the above conditions are satisfied. Thus, the potential 
A^ can be eliminated by the zero of phase for tp~ and it does not represent a real, physical 
electromagnetic field. 

Proceeding in the same way for delivers 

At = -{d^C-^ 
AX = 

M+ = \e\ (40) 

under the conditions d_9 = and 

X'DC = -{d. + d+)^ = -Tdtl (41) 



As above, these conditions makes the electromagnetic tensor vanish. Thus, A'^ does not 
represent a real, physical electromagnetic field. 
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We thus conclude that there are indeed situations for which one or the other autonomous 
Klein-Gordon equation follows from a variational principle, but that these situations are not 
particularly interesting, at least as far as the dynamics of the retained component of \1/ is 
concerned. In the other situations, each Klein-Gordon equation contains terms which induce 
time-irreversibility and the equations do not follow from autonomous variational principles. 
This is not so surprising after all. It only tells us that the two components of \1/ are really 
coupled by the quantum walk and cannot be decoupled without paying the price i.e. loosing 
symmetry and the possibility of writing simple variational principles. 

V. CONCLUSION 

A. Summary 

We have identified the 1-jets of ID discrete-time quantum walks which admit a continuous 
limit. The continuous limit is described by a simple Dirac-like equation which derives from 
a variational principle. Neither the dynamics of the discrete-time walk nor the dynamics of 
its continuous limit are invariant under SU(2) , but the continuous dynamics nevertheless 
belongs to a family of Dirac-like equations which exhibit an elegant symmetry. In the 
continuous limit, each component of the wave function also obeys an autonomous Klein- 
Gordon equation, but these equations do not generally derive from variational principles. 

B. Discussion 

Let us now discuss how these results connect with the literature. As mentioned earlier, 
there are two types of quantum walks, called respectively discrete time^"^ and continuous 
timei'i^'2^ quantum walks. Space is discrete in both types of walks. It is shown in^ that, at 
least formally, discrete time quantum walks become continuous time quantum walks when 
the time step tends to zero. The correspondence between discrete and continuous time 
quantum walks is also explored in^^. Our work investigates what happens when both time 
and length steps tend to zero. The problem was first touched upon by Feynman^, who was 
in search of a simple dynamics which, when quantized, would give back the Dirac equation. 
Feynman's approach is discussed rapidly and the continuous limit of quantum walks 

is dealt with extensively ii^^-^. 
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Reference8^'2S propose continuous limits for several discrete walks with constant coeffi- 
cients, including the Hadamard walk. This seems to contradict the material presented in 
Section 2 above, where we have proved by a rather trivial argument that the Hadamard walk 
cannot have a formal continuous limit. This apparent contradiction is resolved by two re- 
marks. First, the differential equations proposed in^i22 are not obeyed by the original field 
but rather by two composite fields which mix ip' and ip^ in a time dependent manner. Sec- 
ond, the authors of referencea^i^S do not use a strict concept of continuous limit. Indeed, the 
differential equations proposed in^i22 are not the strict formal equivalent of the random walk 
equations when both time- and length steps tend to zero; the differential equations proposed 
in^i22 rather admit discretizations which start with terms identical to those found in the dis- 
crete quantum walk dynamics. This article therefore does not contradict reference^i^', but 
rather complements them by presenting, when it exists, the strict formal continuous limit of 
discrete time-and space-dependent walks in the form of differential equations obeyed by \E'. 

Referencesi^i^^ address only the special case where 9 and C, identically vanish and is a 
constant. These references mention that a Dirac-like equation obeyed by \1' can be formally 
obtained from a discrete quantum walk when the angle ^ tends to zero with the time and 
length steps, but they offer little detailed computations and do not mention the existence 
of Klein-Gordon equations. ReferenceA^ only addresses situations where ^ = C = and 
^ is a constant. This reference does present a detailed computation and both Dirac and 
Klein-Gordon equations are obtained, but the crucial condition that 9 must tend to zero 
together with the time and length step is overlooked. As a result, the equations presented 
iv^ contain various trigonometric functions of 9. Once the correct scaling law is applied, 
these equations turn out to be identical to the power laws and logarithm appearing in our 
work. 

The following four points are not discussed in the previous literature and have thus been 
addressed here for the first time: (i) general situations where none of the the Euler angle 
vanishes and where these angles depend on time and/or space (ii) the existence of several 
scaling laws (iii) the existence of variational principles for the continuous limit (iv) the 
symmetry properties of the continuous limit. The present article is thus a clear on the 
current literature. 

This work can be extended in various directions. The above material makes clear that 
the scaling studied in this article is not the only one which allows for a continuous limit 
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of the discrete walk. Other scahngs, though not as rich as the one retained in this article, 
should naturally be studied in their own right. Detailed numerical simulations should also 
be performed to fully understand the dynamics of the continuous limits. It would also be 
very interesting to investigate the continuous limit of discrete quantum walks defined over a 
higher dimensional space-times and/or quantum walks driven by quantum coins acting on 
higher dimensional Hilbert spaces. And one can only wonders how discrete quantum walks 
can be coupled to space-time geometry and if this coupling can be of practical interest, for 
example in astrophysics or cosmology. It should finally prove very rewarding to investigate 
how decoherent continuous quantum walks connect with relativistic stochastic processes^*^. 
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